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Abstract. Relativistically covariant form of equation of motion for arbitrarily shaped 
dust particle (neutral in charge) under the action of electromagnetic radiation is derived 
\Q , - emission, scattering and absorption of radiation is considered. The result is presented 
in the form of optical quantities used in optics of dust particles. The obtained equation 

o : 

of motion represents a generalization of the Poynting-Robertson (P-R) effect, which is 

O ■ 

"■■■«» standardly used in orbital evolution of dust particles in astrophysics. 
Oh' Simultaneous action of electromagnetic radiation and gravitational fields of the central 
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body - star - on the motion of the particle is discussed. 
Key words: relativity theory, cosmic dust 



1. Introduction 

The first case closely connected with relativistic equation of motion for a particle under 
the action of electromagnetic radiation was presented by Einstein (1905), who calculated 
the change of energy and light pressure for plane mirror. Later on, Robertson (1937) has 
presented equation of motion for perfectly absorbing spherical particle - this case has 
been applied to astronomical situations for several decades as the well-known Poynting- 
Robertson (P-R) effect. 

It is inevitable to have in disposal general equation of motion. Equation of motion for 
all types of particles, not only for the special cases treated by the above cited authors. 
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Real particles scatter radiation in a more complicated manner and particles of various 
optical properties occur. Moreover, thermal emission has to be included, in order to be 
the effect of radiation complete. Thus, this paper generalizes the results presented in 
Klacka (2000a, 2000b, 2001, 2002a) and Klacka and Kocifaj (2001), in which also results 
of detailed numerical calculations for (sub) micron sized dust particles in the Solar System 
may be found - motion of real particles may significantly differ from that corresponding 
to the P-R effect. 

2. Proper reference frame of the particle stationary particle 

The term "stationary particle" will denote particle which does not move in a given iner- 
tial frame of reference. Primed quantities will denote quantities measured in the proper 
reference frame of the particle. 

The flux density of photons scattered into an element of solid angle dft' — sin 9' dd' d<j>' 
is proportional to p'(9',<j>') dCl', where p'{6' \<j>') is "phase function". Phase function de- 
pends on orientation of the particle with respect to the direction of the incident radiation 
and on the particle characteristics; angles 9', </>' correspond to the direction (and orienta- 
tion) of travel of the scattered radiation, 6' is polar angle and it equals zero for the case 
of the travel of the ray in the orientation identical with the unit vector e[ of the incident 
radiation. The phase function fulfills the condition 



The momentum of the incident beam of photons which is lost in the process of its 
interaction with the particle is proportional to the cross section C' ext (extinction). The 
part proportional to C' abs (absorption) is completely lost and the part proportional to 
C' ext — C' abs = C' sca (scattering) is again reemitted. The differential scattering cross 
section dC' sca /dQ' = C' sca p'(9', cf)') depends on the polarization state of the incident light 
as well as on the incidence and scattering directions (e. g., Mishchenko et al. 2002). 

The momentum (per unit time) of the scattered photons into an element of solid angle 




(1) 



dSV i 



is 



dp' sea = - S' C' sca p'{6',4>') K' dSH 



(2) 



where unit vector in the direction of scattering is 



K' = cos 9' e[ + sin 9' cos </>' e' 2 + sin 9' sin <j>' e' 3 . 



(3) 
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S' is the flux density of radiation energy. The system of unit vectors used on the 
right-hand side of the last equation forms an orthogonal basis. The total momentum (per 
unit time) of the scattered photons is 

p' sca = - S'C' sca f p'(6',ct>')K' dfi' . (4) 

c Jin 

The momentum (per unit time) obtained by the particle due to the interaction with 
radiation is 



dt> ' 



S> {c' ext ei - C" sca J K' dfi'j . (5) 



As for the energy we suppose that it is conserved: the energy (per unit time) of the 
incoming radiation, equals to the energy (per unit time) of the outgoing radiation (after 
interaction with the particle). We will use the fact that t' — r, where r is proper time. 

For the sake of brevity we will use dimcnsionlcss efficiency factors Q' x instead of cross 
sections C' x : C' x — Q' x A', where A' is geometrical cross section of a sphere of volume 
equal to the volume of the particle. Equation (5) can be rewritten to the form 



-f- = - S> A' {[Q' ext < cose' > Q' sca ] ei + 



dp' _l 

d t c 

[ - < sin^' cos^' > Q' sca ] e' 2 + [ - < sinfl' sin0' > Q' sca ] e' 3 } . (6) 



Summarizing important equations, we can write them in a short form 

^- = - S' A' {Q[ e[ + Q' 2 e' 2 + Q' 3 e' 3 } ; ^ = . (7) 
arc " (It 



3. Stationary frame of reference 

By the term "stationary frame of reference" we mean a frame of reference in which 
particle moves with a velocity vector v = v(t). The physical quantities measured in the 
stationary frame of reference will be denoted by unprimcd symbols. 

Our aim is to derive equation of motion for the particle in the stationary frame 
of reference. We will use the fact that we know this equation in the proper frame of 
reference - see Eqs. (7). We have to use Lorentz transformation for the purpose of making 
transformation from proper frame of reference to stationary frame of reference. 

If we have a four- vector A^ = (A , A), where A is its time component and A is its 
spatial component, generalized special Lorentz transformation yields 

A ' =7(4° - v-A/c) , 

A' = A + [(7 - 1) v ■ A/v 2 - 7 A°/c] v , (8) 
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where 7 = 1/ \J\ — v 2 /c 2 . The inverse generalized special Lorentz transformation is 
A = 7 (A ' + v ■ A'/c) , 

A = A' + [(7 - 1) v ■ A'/v 2 + 7 A°'/c] v . (9) 
Well-known four-vector is four-momentum: p M = (p°,p) = (E/c,p). 

3.1. Transformation of equation of motion 

Eqs. (9) yield for our case represented by the second of Eqs. (7): 
dE_ _ djS_ 

dT ~ 7 V ' dT ' 

dp dp' ( dp'\ v 



1Z = — + (7 -!)(«• — ) 72 • ( 10 ) 



5. 2. Flita; density of radiation energy 

Let us consider that a beam of incident photons is characterized by flux density of radia- 
tion energy S' = n'hv'c, or S — nhvc, where n' and n are concentrations of photons in the 
corresponding reference frames. The incident photon is characterized by four-momentum 
■p^ = (hv, hue\). Application of Eqs. (8) to the four-momentum of the photon yields: 

V = V W\ , 

e'l = — {ei + [(7 - 1) v-ejv 2 - 7 /c] v} , (11) 

where abbreviation w\ = 7 (1 — v ■ ei/c) is used. Continuity equation c^j^ 1 = 0, where 
= (cn, cne\) is four-vector of current density yields n' — nw\. Inserting this result 
together with the first of Eqs. (11) into the flux density of radiation energy S', we obtain 

S' = wlS . (12) 

3. 3. Transformation of unit vectors 

We have already found transformation of the vector e\. We have to found also transfor- 
mations of vectors e 2 and e 3 . Physics of the corresponding transformation is that the 
vectors e 2 and e 3 describe directions of propagation of light after interaction with the 
particle (see Klacka 2000a for more details). Thus, aberration of light exists for each of 
these unit vectors. Considerations analogous to those for vector ei immediately yield 

e'j = — {ej + [(7 - 1) v-ej/v 2 - j/c] v} , j = 1,2,3, 

Wj 

Wj = 7 (1 - v ■ ej/c) . (13) 
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3.4- Equation of motion 

Inserting Eqs. (7), (12) and (13) into Eqs. (10), one easily obtains 

3 

^ = E Q'jwiSA' (1/ Wj - 7 ) , 



rE^^t-^-^). d4) 

Eqs. (14) may be rewritten in terms of four-vectors: 
dp» wjS A' * 



d T c 

3 = 1 



(15) 



where p M is four- vector of the particle of mass m, = m u M , four- vector of the 
world- velocity of the particle is u M = (7 c, 7 v). We have also other four- vectors 

b» = (l/w 3 ,e 3 /w 3 ) , j -1,2,3. (16) 

It can be easily verified that Eq. (15) yields d m/d r = 0. 
As for practical calculations we introduce: 

b°j = 1M=7(1 + v-e'j/c), .7 = 1,2,3, 

bj = e j / w j = e ^ + 1(7 - 1) w ' e'j I v 2 + 7 / c] v . (17) 
Within the accuracy to the first order in v/c, Eq. (15) yields 



d v S A' ' 



J2 Q'j K 1 - 2vei/c + v-ej/c) e, - v/c] . (18) 



dt m c 

Eq. (15) yields as special cases the situations discussed in Einstein (1905) and Robert- 
son (1937) - Robertson's case is obtained simply putting Q[ — 1, Q' 2 — Q' 3 — 0, Einstein's 
results require a little more calculations. 

3.5. Heuristic derivation 

Since we completely understand physics of Eq. (15), we are able to present short simple 
derivation of Eq. (18), now. We have: 

i) dv'/dt = [S'A'/(mc)\ £-=i Q[ e[ (see Eq. (7)); 

ii) S' = 5(1 — 2v ■ ei/c) (see Eq. (12)), due to the change of concentration of photons 
n' = n(l — v ■ ei/c) and Doppler effect v' = v{\ — v ■ ei/c); 

hi) e'j = (1 + v ■ ej/c)ej — v/c, j = 1, 2, 3 (aberration of light - see Eq. (13)). 
Taking into account these physical phenomena, we finally obtain Eq. (18). 
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4. Thermal emission 

If the particle's absolute temperature is above zero, it can emit as well as scatter and 
absorb electromagnetic radiation. 

4-1. Proper reference frame of the particle - stationary particle 

This subsection is based on the presentation by Mishchenko et. al. (2002). 

The emission component of the radiation force acting on the particle of absolute 
temperature T' is 

F' e (T') = -- r du' [ r' K' e (r',T>,J) dr' , (19) 

c JO Ji-n 

where the unit vector r = r 1 /r' is given by position vector r' of the observation 
point with origin inside the particle (the emitted radiation in the far-field zone of the 
particle propagates in the radial direction, i. e., along the unit vector r'), ui' is (angular) 
frequency of radiation, 

K' e (r',T',J)=I' b {T',J) {k' u (f',u/) - j[ Z'^{r',r",J) df"} , (20) 
and the Planck blackbody energy distribution is given by the well-known relation 

The functions K' n and Z' n 
K' n = |f Jm(S^ + S' 22 ) , 

Z'n = (\S'n\ 2 + \S[ 2 \ 2 + \S' 21 \ 2 + \S' 22 \ 2 ) / 2, (22) 
are given by elements of the amplitude matrix S'(n sca , n mc ): 

E lsca {r' h' sca ) = ^-J— - S\ti sca , h' mc ) E'™ c , r' 00 . (23) 

We assume that the incident field is a plane electromagnetic wave given by 

E' mc (r') = E'j nc cxp (i fc 1 n ,mc • r') , J 2 = c 2 \k\\ 2 , (24) 

since we consider vacuum as a homogeneous nonabsorbing medium. We use a couplet 
( , ip' ), where G ( 0,7r ) is the polar (zenith) angle measured from the positive 
z'— axis and ip' G ( 0, 2 tt ) is the azimuth angle measured from the positive a;'— axis in 
the clockwise direction when looking in the direction of the positive z 1 — axis. Since the 
medium is assumed to be nonabsorbing, E' = E'#, 1? + E' , ip' , where E'$, ■& lies in the 
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meridional plane (i. e., the plane through n and the z' — axis), whereas the component 
E' , (p is perpendicular to this plane; & and Cp are the corresponding unit vectors such 

a / nf a / 

that n = i9 x ip : 

E'™ c = E'^, c ti mC + E'™f Cp' mc , 

E ,sca = gsca ft*™ + gsca £sca ^ ^ . ^sca^,^ _ q ^ / _ qq) . ( 2 5) 

^.2. Stationary frame of reference 

Using orthonormal vectors e[, e' 2 and e' 3 defined by Eq. (3), we can write Eq. (19) in the 
form 

F' e (T') = £ [F' e {T')-e' t ] ej. (26) 
Using four-vectors defined by Eq. (16), we can immediately write 



d 

where 



E Q'eJ V > (27) 



Qej = -e' r - [ dJ [ f ' K' e (f , T', u/) df , j = 1,2,3. 

c JO J4tt 



As a consequence, Eqs. (16) and(27) yield 



i=i 

As a consequence, dm/ d r = 0. 



(28) 



77 j =E^/ C - ( 29 ) 
Mass of the particle decreases due to the thermal emission, alone. 



5. Electromagnetic interaction total 

5.1. Proper reference frame of the particle - stationary particle 
We can write, on the basis of Eqs. (7), (26) and (28): 

j = l 

5. 2. Stationary frame of reference 

On the basis of Eqs. (12), (16) and (30) we obtain 



(31) 
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6. Gravitational and electromagnetic radiation fields 

Let us consider orbital evolution of real dust particle in the Solar System, under ac- 
tion of gravitational and electromagnetic radiation fields of the Sun. We can write (also 
formulation presented in Klacka 2002b is used) 



dv = _ GMv + GMv » r, _ ^ v _e^ + 

at r z r z ^-^ W c c / 



where 



E Qei [(l + V ~^) * - > ( 32 ) 



ej = (1 — v-e'j/c)e'j + v/c, j — 1,2,3, 



/3i = 



© 



G M m c 



/ B (A) {C' ext {X/w) C' sca (X/w) g[(X/w)} dX , 
Jo 



f3 2 = „ ' A ; / S (A) { - C' sca (X/w) g' 2 (X/w)} dX , 

= CM Zir / ^ - C aca(A/^ ^(V™)} dA , 

(_r Mq m c j 

w = I — v ■ e\ / c , (33) 

Rq denotes the radius of the Sun and -B©(A) is the solar radiance at a wavelength of A; 
G, Mq, and r are the gravitational constant, the mass of the Sun, and the distance of the 
particle from the center of the Sun, respectively. The asymmetry parameter vector g' is 
defined by g' = {1/C' sca ) J n' (dC' sca / dQ')dfl' , where n' is a unit vector in the direction 
of scattering; g' = g[ e[ + g' 2 e' 2 + g' 3 e' 3 , e[ = (1 + v ■ e\/c) e\ - v/c, e- • e'j = Sij. 

7. Conclusion 

We have derived relativistically covariant equation of motion for dust particle under 
the action of electromagnetic radiation - see Eq. (31). As a generalization requires, 
Eq. (31) may be reduced to special cases treated by Einstein (1905), Robertson (1937) 
and Klacka (2000a). Application of Eq. (31) to Solar System is represented by Eqs. 
(32) and (33). Some other accelerations can be added to the right-hand side of Eq. 
(32) - e. g., gravitational perturbations of planets, solar wind effect [(dv /dt) sw = 
— (r]f3iGMQ/(Q' 1 r 2 ))((v ■ e 1 /c)e 1 + v/c);t] « 0.3], or some other nongravitational ac- 
celerations. 
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